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Abstract 

The main goal of this paper is to get in a straightforward form the field equations in metric f(R) 
gravity, using elementary variational principles and adding a boundary term in the action, instead 
of the usual treatment in an equivalent scalar-tensor approach. We start with a brief review of the 
Einstein-Hilbert action, together with the Gibbons- York-Hawking boundary term, which is mentioned 
in some literature, but is generally missing. Next we present in detail the field equations in metric 
f(R) gravity, including the discussion about boundaries, and we compare with the Gibbons- York- 
Hawking term in General Relativity. We notice that this boundary term is necessary in order to have 
a well defined extremal action principle under metric variation. 
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1 Introduction 

General Relativity (GR) is the most widely accepted gravity theory proposed by Einstein in 1916, and 
it has been tested in several field strength regimes being on of the most successful and accurate the- 
ories in physics [1]. The field equations can be obtained using a variational principle, from the well 
know Einstein-Hilbert action [2]- [6]. The methodology leads to a boundary contribution which is usually 
dropped out [7], [8], setting null fluxes through Gauss-Stokes theorem. It can be done by imposing that 
the variation of the metric and its first derivative vanishes in the boundary [3]. These conditions can 
be relaxed whether a boundary term is introduced, called the Gibbons- York-Hawking boundary term 
[9], [10]. With this boundary term is necessary only to fix the variation of the metric in the boundary. 
There are some references [3]- [5] where this boundary term is shown explicitly. 

However, GR is not the only relativistic theory of gravity. In the last decades several generalizations 
of Einstein field equations have been proposed [11]- [14]. Within these extended theories of gravity nowa- 
days a subclass, known as f(R) theories, are an alternative for classical problems, as the accelerated 
expansion of the universe, instead of Dark Energy and Quintessence models [15]-[21]. f{R) theories of 
gravity are basically extensions of the Einstein-Hilbert action with an arbitrary function of the Ricci 
scalar R [22]- [25]. The field equations were founded in [26], and including boundary terms in fourth order 
gravity in [27]-[29]. The Gibbons- York-Hawking like term in f{R) gravity was explored in [30], with an 
augmented variational principle in [31], [32], and using a scalar-tensor framework in [33]- [38]. Here we ob- 
tain the field equations from a metric f(R) action with boundary terms, using only variational principles. 
We get a well constrained mathematical problem setting Sg a p = and SR = in the boundary. 
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2 General Relativity: The Einstein-Hilbert action with the Gibbons- 
York-Hawking boundary term 

We consider the space-time as a pair (A4,g) with M. a four-dimensional manifold and g a metric on 
A4 . GR is based on the Einstein's Field equations (without cosmological constant and geometrical units 
c = 1), which gives the form of the metric g a p on the manifold M.: 

Ra/3 — yRgaf! = K>T a p, (2.1) 

where R a p = R^p is the Ricci tensor, R = R a @R a p the Ricci scalar, and T a p the stress-energy tensor, 
with k = 87rG, and sign convention (— , +, +, +). The Riemann tensor is given by: 

n P-yS — a 7 i 5/3 — °S L 7 /3 + 1 7a 1 5/3 1 aS l j0, \ Z - A > 

in terms of the connections Tp* . The Einstein field equations can be recovered by using the variational 
principle SS = 0, with S expressing the total action. In terms of Einstein-Hilbert action Seh, Gibbons- 
York-Hawking boundary term Sqyh and the action associated with all the matter fields Sm, the total 
action can be written by [4]: 

S = —(Seh + Sgyh) + Sm, (2-3) 

where 

Seh = ! d 4 x^g~R, (2.4) 



J v 

S GYH = 2l d 3 yej\h\K, (2.5) 

JdV 

here V is a hypervolume on VW, dV its boundary, h the determinant of the induced metric, K is the trace 
of the extrinsic curvature of the boundary dV, and e is equal to +1 if dV is timelike and —1 if dV is 
spacelike (it is assumed that dV is nowhere null) . Coordinates x a are used for the finite region V and y a 
for the boundary dV. Now we will obtain the Einstein field equations varying the action with respect to 
g a P '. We fixed the variation with the condition [3], [4] 



J a 13 



= 0, (2.6) 

dV 



i.e., the variation of the metric tensor vanishes in the boundary dV . We use the results [4], [7] 

Sgafi = -g^gpJtT, 5g a0 = -sT^Sg^, (2.7) 

SV^g = -\v^gg a pS g ^, (2.8) 

SR^ s = ^(6TJ p )-V s (6T^p), (2.9) 

SR a p = V 1 (ST} a )-Vp(ST'r a ). (2.10) 

We give a detailed review for the variation principles in GR following [3], [4] and [7],. The variation of 
the Einstein-Hilbert term gives 

SSeh = I d i x(RS^g + ^gSR). (2.11) 
Jv 

Now with R = g al3 R a f3, we have that the variation of the Ricci scalar is 

SR = 5g a0 R a p + g ap 5R a p. (2.12) 
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using the Palatini's identity (2.10) we can write [7]: 

SR = 5g a ' 3 R a0 + 9^(V^6T} a ) - V^I^)), 
= 5g^R a0 + V a (g al3 (Sr% a ) - g^ST*)), 



(2.13) 



where we have used the metric compatibility Vjg a p = and relabeled some dummy indices. Inserting 
this results for the variations in expression (2.11) we have: 

SS BH = [ d 4 x(R5 v ^g + v^gSR), 
Jv 

= J v d 4 x (—Rg^V^gSg^ + RapV^Sg^ + V^V ^ (ST^ a ) ~ g^iSTZ^ , 

= - \Rg a ^jSg^ + j ^x^gsf a {g a ^sr% a ) - g^iST^)). (2.u) 



Denoting the divergence term with SSb, 

ss b = f d'^v^^^j-r^)), 

Jv 

we define 

V- = g^(ST^ a )-g^(STl 7 ), 



(2.15) 
(2.16) 

(2.17) 



then the boundary term can be written as 

SS B = I d 4 xJ~gV a V a . 
Jv 

Using Gauss-Stokes theorem [4], [7]: 

f d n x^\g~\V^ = I <r- 1 yey/\h\n ll A>>, (2.18) 

JV JdV 

where is the unit normal to dV. Using this we can write (2.17) in the following boundary term 

SS B = I d 3 ye^/\h\n a V\ (2.19) 

JdV 

with V a given in (2.16). The variation STp a is obtained by using that Tp a is the Christoffel symbol {^ a }' 



(2.20) 



getting 



^ r /3a = s Qff CT7 [ d P9ya + d a g~«f3 - d 1 g f3a ] j . 



= \ 6 V ai [ d P9~fa + dagyfs - d 1 gp a ] + \^g ai [dp(Sg ja ) + d a {Sg lj3 ) - d^{5gp a )] . 

From the boundary conditions Sg a p = Sg a * = the variation (2.21) gives: 

1 



ST 



/3a 



dV 2 



g" 1 [dp(5g ia ) + d a (Sg 7 p) - d 1 (5g ! j a )\, 



and 



n a/3 



i)V 



-g M [dp{5g^ a ) + d a {5g lf} ) - d 7 (Sg Pa )] 



-g vl d a {Sg ul ) 



(2.21) 

(2.22) 
(2.23) 



wc can write 



Vn 



dv 



9 '<j p. 9 



ad 



-g^ 1 [dp{8g ia ) + d a (Sg lfl ) - d 7 (6gp a )] 



gcr^g 



-g vl d a {8g vl ) 



]f a g^ [8p(8g ja ) + d a (5g lf} ) - £,(%»)] - \^g v "> [d a (Sg VJ )] , 



= g aP [dp(5g„ a )-d a (5gp a )]. 
We now evaluate the term n a V a \ dv by using for this that 



then 



n°V n 



i)V 



g a P = h ap +en a n p , 

n°{h a(i + en a n p )[dp[8 g< r a ) - d^gp*)], 
= n' J h a0 [dp{Sg lja )-d <J {Sgp a )}, 



(2.24) 
(2.25) 

(2.26) 



where we use the antisymmetric part of en a rfi with e = n^rin = ±1. To the fact 5g a p = in the 
boundary we have h a " dp(8g aa ) = [4]. Finally we get 



= -rfh^d^g^). 



Thus the variation of the Einstcin-Hilbcrt term is: 



8S EH = J ^xV^g^p - ^Rg a ^Sg al3 - j d 3 ye^\h\h al3 d a (Sgp Q )n 
Now we consider the variation of the Gibbons- York-Hawking boundary term: 

SSgyh = 2 f (Pye^/\h\5K. 

JdV 



(2.27) 



(2.28) 



(2.29) 



IdV 

Using the definition of the trace of extrinsic curvature [4]: 

K = V a n a , 

= {h af} + en a n /3 )V p n a 

= h^V fj n a , 

= h afi (di,n a -ri a r h ), 



(2.30) 



the variation is 



SK = -h afi 5r 



= -\h a ^g^[d fj {5g. a ) + d a {5g aP ) - d a {5gp a )]i 
= - \^ a(i [9/3 (Sgaa ) + d a (Sg al3 ) - d a (Sgpa )] rf , 
= \h^d„{8g Pa )n°. 



(2.31) 



This comes from the variation ST^ a evaluated in the boundary, and the fact that h a P dp{8g aa ) = 0, 
h a P d a (8g a p) = 0. Then we have for the variation of the Gibbons- York-Hawking boundary term: 



SS, 



GYH 



dv 



d 3 ys^/\h\h af) d G {8gp a )n° 



(2.32) 
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We see that this term exactly cancel the boundary contribution of the Einstein-Hilbert term. Now, if we 
have a matter action defined by: 



Sm 



/ (Px^f=gCM\goci},i>\, 
Jv 



where tp denotes the matter fields. The variation of this action takes the form: 

SS M = / d 4 xS(y/^g£ M )i 
Jv 



= / d A x \f—g 
Jv 

as usual, defining the stress-energy tensor by: 

n dC M , 



dC 



dg a P 2 



„o>/8 



T a p 



'dg 



a/3 



2 SS 



M 



-9 < 



then: 



SS M = -- / d 4 xy/^gT a pt 



v 



(2.33) 



(2.34) 



(2.35) 



(2.36) 



imposing the total variations remains invariant with respect to 5g a ' . Finally the equations are writing 
as: 

1 SS 1 

; — T-Za = 0, R a p - 7,R9ap = KT a p, (2.37) 
y/-g dg a P 2 

which corresponds to Einstein field equations in geometric units c = 1. 



3 Field equations in f(R) gravity 

As we mentioned above the modified theories of gravity have been studied in order to explain among the 
accelerated expansion of the universe. One of these theories is the modified f(R) gravity which consists 
in add additional higher order terms of the Ricci scalar in the Einstein-Hilbert action [15], [16], [25]. There 
are three versions of f(R) gravity: Metric formalism, Palatini formalism and metric-affine formalism [25]. 
Here we focus only in the metric formalism; for a detailed deduction of field equations in the Palatini 
and the metric-affine formalism see [39], [40]. Again, we consider the space-time as a pair (J\4,g) with 
M. a four-dimensional manifold and g a p a metric on M . Now the lagrangian is an arbitrary function of 
the Ricci scalar £[g a ^] = f(R), the relation of the Ricci scalar and the metric tensor is given assuming a 
Levi-Civita connection of the manifold, i.e. a Christoffel symbol. This lagrangian was presented in [31] 
using augmented variational principles. The general action can be written as [38]: 

Smod = (Smet + S'gY h) + ^Mt (3-1) 

with the bulk term 

S met = [ d 4 x^f(R), (3.2) 
Jv 

and the Gibbons- York-Hawking like boundary term [30], [38] 

S GYH =2<f d 3 yey^\h\f'(R)K : (3.3) 

JdV 
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with f'(R) = df(R)/dR. Again, Sm represents the action associated with all the matter fields (2.33). 
We fixed the variation to the condition 



lap 



= 0. 



(3.4) 



HV 



(3.5) 



(3.6) 



First, the variation of the bulk term is: 

SS met = [ d 4 x (f(R)5y/=g + y/=g6f{R)), 
Jv 

and the functional derivative of the f(R) term can be written as 

Sf(R) = f'(R)SR. 

Using the expression for the variation of the Ricci scalar: 

SR = SgTPRcp + S7 a (g^(ST^ a ) - g a °(6TZ 7 )), (3.7) 

where the variation of the term g OL P(8T'p a ) — g aa (8T1 tl ) is given in A. With this result the variation of 
the Ricci scalar becomes 

SR = Sg^R aP + Va(g a P(6T% a ) - ^(ST^)), 
= 5g af} R al3 + g^V a V a {5g^) - V CT V 7 (^), 

= SsT^Rafi + g*pO(5g a P) - VaVpiSg"?). (3.8) 

Here we define □ = V a V a and relabeled some indices. Putting the previous results together in the 
variation of the modified action (3.5): 



8S n 



d 4 x {f{R)5^~g + V~g f'(R)SR) , 



= j^d 4 x (-.f(R)\^g~g aP 5g afi + f ' {R)y/=g{Sg^ R a p + g a pU{5g a ^) - V a V p (5g afi ) 

= J^d 4 x^(f'(R)(Sg^R a p+g a pD(5g^) - V„V/j(5^)) - f{R)\g a p8g^ 

Now we will consider the next integrals: 

/ #xy/=gf(R)g a pU(8tf*), [ d*x V=fl/'(iJ)V„V ?(5g a0 ). 
Jv Jv 

We shall see that these integrals can be expressed differently performing integration by parts. For this 
we define the next quantities: 



(3.9) 



(3.10) 



M T = f(R)g a ^ T (Sg a0 ) - 5g^g a pW T (f'(R)), 



and 



(3.11) 
(3.12) 



iV CT = /'(i?)V 7 (^) - Sg^V y (f(R)). 

The combination g aT M T + N" is 

sTMr + N° = f'(R)g a ^(Sg^) - 8g^g a ^(f'(R)) + f'(R)V 7 (8g^) - <5.g^V 7 (/'(i?)), (3.13) 

in the particular case f(R) = R, the previous combination reduces to the expression (2.16) with equation 
(A. 11). The quantities M T and N" allow us to write the variation of the bulk term (3.9) in the following 
way (for details see B): 

SS met = J^d i x^^f'(R)R a p+g af! af'(R)-\7 a \7 p f(R) - f{R)^g a ^5g afi 

+ I d 3 y£^J\h\n T M T + I d 3 y £^/\h\n r7 N a . (3.14) 

JdV JdV 
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In the next section we will work out with the boundary contribution from (3.14), and show how this 
terms cancel with the variations of the S' GY h action. 

3.1 Boundary terms in f(R) gravity 

We express the quantities M T and N a calculated in the boundary dV. Is convenient to express them in 
function of the variations 8g a j3- Using the equation (2.7) in (3.11) and (3.12) yields : 



and 



M T = ~f'(R)g af3 V T (5g aP ) + g^Sg a ^ T (f'(R)), 
N° = -f{R)g^g^ 1 {8g IJlv ) + . 9 CT V l/ %,V 7 (/'(i?)). 



(3.15) 



(3.16) 

To evaluate this quantities in the boundary we use the fact that 6g a p\dv = 8g a ^\ev = 0, then the only 
terms not vanishing are the derivatives of 8g a p in the covariant derivatives. Hence we have 



■f'(R)9 af, d r (6g aP ), 



and 



(3.17) 



(3.18) 



We now compute n T M T \„ v and n a N a \„ v which are the terms in the boundary integrals (3.14) 



n T M T 



= -f'(R)n T (en a nP + h a ?)d T (Sg a p), 



= -f{R)rfh a ^d a {5g a p) 
where we rename the dummy index r. In the other hand 



(3.19) 



ov 



= -f'{R)n a {h' 7 ^ + erfn^K 1 " + erfn^d^Sg^), 

= -f(R)n^(h^ + erfn v )d^Sg^), 
= -f(R)n»h^d 7 (8g^) 
= 0, 



(3.20) 



where we have used that n a h atl = 0, e 1 = 1 and the fact that de tangential derivative h lv d 1 (8g l _ LU ) 
vanishes. With this results the variation of the action S me t becomes: 

8S met = J^x^( K f l {R)R a p+g a pUf'(R)-V a Vpf'{R) - f(R)^g a p)8g af3 

- I d 3 ye^\h\f{R)n°h af) d a {8g afj ). (3.21) 

JdV 

We proceed with the boundary term S' GYH in the total action. The variation of this term gives 



6S' r 



GYH 



d 3 yeV\h\(Sf'(R)K + f'(R)8K), 



2 <f d 3 ye v / \h\(f"(R)8RK + f{R)8K) 

JdV 



(3.22) 
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Using the expression for the variation of K 1 equation (2.31), we can write 



5S' GYH = 2£^d 3 yeVW\^r(R)SRK + ^f(R)h^d a (6 gf]a )n^, 

= 2/ d 3 ye^/\h\f"(R)SRK + <f d 3 y e^\h\f'(R)h af3 d a (5g Pa )n° . (3.23) 
Jav Jav 



We see that the second term in (3.23) cancels the boundary term in the variation (3.21), and in addition 
we need to impose 5R — in the boundary. Similar argument is given in [38]. 

Finally, with the variation of the matter action, given in (2.36), the total variation of the action of 
modified f(R) gravity is: 



SS mod = ^Jd A xyf^(f{R)R a p+g a0 Uf(R) - V Q V (3 /'(i?) - ^f(R) g Q ^Sg ap 

-ljd 4 x^T a0 Sg^. (3.24) 



Imposing that this variation becomes stationary we have: 



1 SS mo d 



1 



^ = f(R)R a p+g af) nf'(R) - V a V p f\R) - -f(R)g a p = «T a/3 , 



(3.25) 



which are the field equations in the metric formalism of f(R) gravity. 

4 Conclusions 

We have obtained the field equations in the metric formalism of f(R) gravity by using the direct results 
from variational principles. The modified action in the metric formalism of f(R) gravity plus a Gibbons- 
York-Hawking like boundary term must be written as: 



<? - 1 
2n 



f=g~(f{R)+2KL M [9 a p,i>])+2 I ePyey/\h\f'(R)K 
K 7 Jav 



(4.1) 



with f'(R) = df(R)/dR and Cm the lagrangian associated with all the matter fields. From the quantities 
M a and N a , defined in (3.11) and (3.12) respectively, we recovered GR plus Gibbons- York-Hawking 
boundary term in the particular case f(R) = R- We see that including the boundary term, we have a 
well behaved mathematical problem setting both, 8g a p = and SR = in dV. 
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A Evaluation of the term g a ^{5T a [ja ) - g aa {5Tl 1 ) 

We already have calculated the variation STg a : 

ST 0a = [dpgya + d a g lf3 - djgpa] + ^g aj [d {8g^ a ) + d a {5g lP ) - d 1 {5g Pa )] , (A.l) 

writing the partial derivatives for the metric variations with the expression for the covariant derivative: 

V 7 <5g Q ^ = d 7 5g a p - T° a 8g al) - T^Sg aa> (A. 2) 

and also using that we are working in a torsion- free manifold i.e., the symmetry in the Christoffel symbol 
= T^, we can write: 

5T °3 a = \ 5 9 ai [dpg-ya + d a g l/3 ~ d 1 gp a ] + -g a ~< \V p{5g ia ) + V a (8g lP ) - V 7 (%„) + T^Jg lX + Y X p8g Xl ] , 
= \ 5 g ai [dpg ia + dag^p - d 1 gp a ] + g' 71 T x ia 5g lX + \g ai [V 'p(Sg ia ) + ^a(Sg 7 p) - V 7 (5ff /8a )] , 

(A.3) 

using equation (2.7) in the second term: 

STp a = \jg ai [dpg ya + d a g lP ~ d l90a ] - 5g^g^g^g x „T x a + X -g^ [W p{8g ia ) + V a (5g ll3 ) - V 7 (% a )] , 
= 6g° v g Xv T} a - 5g^5°g x „T x 0a + \g^ [Vp(6g ia ) + V a (Sg 10 ) - V 7 (% a )] , 

= Sg^g Xv T x a - 5g™g Xv T x a + \g^ [V ' p(8g ia ) + Xt a {Sg lP ) - V 7 (% a )] . (A.4) 
Then we have 

s^ 3a = \g"~<[Vf S {5g ai ) + V a (% 7 ) - V 7 (% a )], (A.5) 

and similarly 

6Tl 7 = ^[V a (Sg ai )]. (A.6) 

However it is convenient to express the previous result in function of the variations Sg a @ , we again use 
(2.7): 

= [Vpi-ga^Jgn + Vai-gp^uSgn - V- ( (-g^g a J g n] , 
= ~a^[9a^^p(Sg^)+g^ v ^ a (Sg^)-g p ^g av ^(8g^)], 
= -\[6l9cvV > W") + 5Zg fj ^ a {8gn - g Pfl g a „g^V 7 (6 9 n] , 
= ~[gcn^p{S g r) + gfi^ a (S S r)-gp ll g a uV a {Sgi a ')], (A.7) 
where we write V 7 = g aj V-y. In a similar way: 

^ 7 = -^V a (^)- (A.8) 
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Now we compute the term g aP (ST^ a ) - g Q<T (5r^ 7 ) 

9 aP iST% Q ) - 9 aa (WZ-y) = - l^garrVpitgr) +g aP 9MV a (6g^) - g^g ai N° (Sg^)] (A.9) 

- [g aa g^ a {6 9 n]), 

= - \([%Vf>{8<r)+s!;v a {s<r) - a^vwo] - [v.9 M v a (^)]), 

'[V 7 (5<T) + V 7 (<5 5 ^) - fl ^V'(5«r)] - [<?^VW)] 



then we have, 



2 

= - 1 (2V 7 (<5 5 ^) - 2g^{5 9 n) , (A.10) 
g^(<5r^) - 5 Qff (<5r2; 7 ) = <?^V W) - V 7 (5<T ). (All) 



B Integrals with M r and TV 7 

Taking the covariant derivative in M CT : 

V^A/ T = V r (f{R)g af )VASg af> )) - V T {6g^g a pV T (f'(R))), 

= \7-(fi(R))g ap \7 T (Sg a P) + f{R)g a pU{S^) ~ ^7 T (Sg^)g^V T (f (R)) - 5g^g a0 D(f(R)), 
= .fi(R)g a[i U(5g^) - 5g^g a pU{f{R)). (B.l) 

Here we have used the metric compatibility V T g a /3 = 0, integrating this expression 

/ d A x^g~V T M T = f d i x^df'(R)g aP a(Sg af3 )- f d A x ^gSg^g a0 n(f'(R)), (B.2) 
Jv iv Jv 

using again the Gauss-Stokes theorem (2.18), the first integral can be written as a boundary term: 

/ d 3 ye^/\h\n T M T = f d 4 x (R)g a0 a(5g^)) - f d 4 x ^Sg 01 ? g a/3 n(f'(R)), (B.3) 

JdV JV JV 

then we can write: 

d 4 x^f'(R)g a ^a(Sg afj )= [ d A x V=gSg a/3 g a pD(f (R)) + I d 3 y e^/\h\n T M T . (B.4) 

V JV JdV 

In a similar way. taking the covariant derivative of N a : 

V CT iV CT = V CT (/'(i?)V 7 (<5 5 ^)) - V CT (^V 7 (/'(i?))), 

= V CT (/'(i?))V 7 (<^) + /'(i?)V CT V 7 (^) - V CT (^)V 7 (/'(i?)) - SgTVrWiR)), 

= f{R)V a V p {8g^) - Sg^VrVpif'iR)), (B.5) 

integrating: 

/ d 4 x^V a N° = [ df-xy/=gf\R)V a VfiW)- [ d^x^—gSg^VvVpif'iR)), (B.6) 
Jv Jv Jv 

using again the Gauss-Stokes theorem we can write: 

f d A x^—gf'{R)V a Vp(5g^) = f d*x y/—g~Sg al3 V( r Vp(f'(R)) + I d 3 y e^\h\n a N° . (B.7) 

Jv Jv JdV 
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